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Abstract: In Central Europe, a large portion of post-mining sites were 
afforested with Scots pine, which is characterized by good adaptability 
and a tolerance for poor habitat at the beginning of forest ecosystem 
development. Conversion of monoculture on mine sites into more biodi- 
verse mixed hardwood forests, especially on more fertile deposits, can be 
an emerging need in this part of Europe in next decades. The ability to 
classify the forests at these post-mining sites will facilitate proper species 
selection as well as the management and formation of the developed 
ecosystem’s stability. This work describes the guidelines that can be 
followed to assess reclaimed mine soil (RMS) quality, using the mine 
soil quality index (MSQI) and a classification of developed forest sites as 
a basis of tree-stand species selection and conversion of pine monocul¬ 
tures. The research was conducted on four post-mining facilities (lignite, 
hard coal, sulphur, and sand pit mining areas) on different RMS sub¬ 
strates dominant in Central Europe. Soil quality assessment takes into 
account the following features of the soil: texture soil nutrients (Ca, Mg, 
K, Na, P); acidity (pH KC1); and C or g-to-N, ratio in the initial organic 
horizon. An analysis was conducted of classification systems using the 
MSQI validation correlation (at p =0.05) with vegetation features af¬ 
fected by succession: aboveground biomass of forest floor and ecological 
indicators of vascular plants (calculated on the basis of Ellenberg’s (2009) 
system). Eventually, in the analysed data set, the MSQI ranged from 
0.270 for soils on quaternary sands to 0.720 for a mix of quaternary 
loamy sands with neogene clays. Potential forest habitat types and the 
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role of the pine in the next generation of tree stands on different RMS 
parent rock substrate were proposed. 

Keywords: mining sites, afforestation, pine monoculture conversion, soil 
quality index, forest habitat classification. 

Introduction 

Different components of mining activities including exploration, 
extraction, and processing impose extensive physical, chemical, 
and biological changes on the environment due to the nature and 
characteristics of the activities (Daniels and Stewart 2000; Hiittl 
and Weber 2001; Machaina 2001; Pietrzykowski and Krzak- 
lewski 2007). 

Correct identification of habitat conditions at mining sites and 
planning species composition in afforestation both have a fun¬ 
damental impact on the stability of reclaimed stands and devel¬ 
oped forest ecosystems (Gale et al. 1991; Heinsdorf 1996; Bur¬ 
ger and Kelting 1999; Knoche at al. 2002; Krzaklewski and 
Pietrzykowski 2007). Diagnosis and classification of habitats on 
mining soils reclaimed for forestry is of paramount importance 
for a proper selection of species composition in afforested sites. 
Proper selection of tree species means their ecological require¬ 
ments are appropriately adjusted to habitat conditions (climate 
conditions, fertility and soil humidity). 

In Central Europe, a huge area of reclaimed post-mining sites, 
especially the Lusatia Mine District in Germany and in Central 
Poland, were afforested with a pine monoculture (Heinsdorf 
1996; Knoche 2005; Baumann et al. 2006; Pietrzykowski 2010). 
This forest management and afforestation practice is based on 
the assumption that the initial habitats should follow primary 
succession and, during the first stage, the biotope is colonized by 
pioneers. 

The Scots pine (Pinus sylvestris L.) is a very useful for recla¬ 
mation and afforestation on post-mining sites due to its adapta¬ 
bility and tolerance of poor habitat (Knoche 2005; Baumann et al. 
2006; Pietrzykowski and Socha 2011). This species is native to 
Europe and Asia, ranging from Great Britain and Portugal in the 
west, Eastern Siberia in the east, the Caucasus Mountains in the 
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south and as far north as well inside the Arctic Circle in Scandi¬ 
navia (Farjon 2005). In the case of post-mining sites with more 
fertile deposits, there is a movement toward converting single¬ 
layered Scots pine monocultures into mixed hardwood forests for 
increasing biodiversity, tree-stands resistant to demage caused by 
insects and fungi gradations and better habitats utility for hard¬ 
wood species. Generally, the transformation of Scots pine mono¬ 
culture on potentially fertile soils into close-to-nature mixed 
hardwood forests stands with deciduous species is an important 
forestry practice in Europe in last decades (Buczko et al. 2002). 
This practice accelerates better equilibrium between litter pro¬ 
duction and decomposition and eventually the distribution of 
humus stock in the organic layers and the mineral horizon 
(Fischer et al. 2002). 

In the case of developed ecosystems on post-mining sites, the 
most important issue is their biodiversity and stability, which are 
characteristic of mixed forests. Habitat classification is the basis 
for species selection and stand composition in forestry. In the 
case of post-mining sites afforested with pioneer species like the 
Scots pine, the soil quality index would be a good tool for plan¬ 
ning monoculture conversion to mixed forests in the next genera¬ 
tion of tree stands. 

When planning the reconstruction of species composition the 
best time to start is at the first generation stand, when the as¬ 
sessment of habitat conditions forming under dynamic ecosystem 
succession is possible (Knoche 2005). Carrying out soil quality 
assessments with regard to tree species selection before estab¬ 
lishing the first forest stands on mining sites will also pay off in 
the long run (Pietrzykowski 2010). 

Methods of habitat productivity assessment, developed for 
“natural” forests are of little use in the classification of habitats 
on mining soils reclaimed for forestry (Pietrzykowski and Socha 
2011). Mining and tipping result in the mixing of deposits differ¬ 
ing in (Quaternary, Neogene, Carboniferous strata). This affects 
in diversity of soil texture (mixing sand, clay, silt), and in many 
other parameters (pH, soil fertility and nutrient status, trace ele¬ 
ments, air-water properties) compare to “natural” soils. Thus the 
newly fonned habitats are completely different from “natural” 
habitats. 

The newly-formed soils are characterized by highly change¬ 
able chemical and physical properties and the consequential large 
spatial variability in habitat conditions. Characteristic features of 
reclaimed mine soils (RMS) include: lack of soil organic matter 
(SOM), nutrient deficiency (mainly of nitrogen and phosphorus); 
low pH-values due to acid mine drainage; unfavourable air-water 
properties; and salinity (Daniels at al. 1992; Burger at al. 1994; 
Andrews et al. 1998; Heinsdorf 1996; Katzur and Haubold-Rosar 
1996; Daniels and Stewart 2000). However, elements’ deficiency 
in RMS is not the only problem as there may also be excessive 
concentrations of sulphur from pyrite oxidation and acid mine 
drainage on the Neogene strata (Katzur and Haubold-Rosar 
1996). These features restrict harmonious nutritive conditions for 
trees in recreated forest ecosystems (Burger et al. 1994; Heins¬ 
dorf 1996; Knoche 2005; Knoche et al. 2002, Pietrzykowski 
2008). 

An objective comparison of soil quality on the basis of quanti- 
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tative indices (e.g., soil quality index SQI) is important for sus¬ 
tainable forest management and stability of the ecosystem, which 
also has protective and non-productive functions (connected with 
aesthetics and well-being). A unifonn system of soil quality as¬ 
sessment based on the SQI is also important for monitoring the 
anthropogenic impact on the environment (Schoenholtz et al. 
2000) and is useful for universal assessment of potential forest 
soil productivity (Gale et al. 1991; Burger and Kelting 1999; 
Rnoepp et al. 2000, Brozek et al. 2011). In practice, soil quality 
can be evaluated on the basis of its chemical, physical, and bio¬ 
logical properties (Doran and Parkin 1996; Henderson 1995; 
Regenold and Palmer 1995; Harris et al. 1996; Nambiar 1997; 
Powers et al., 1998; Schoenholtz et al., 2000). 

An example of SQI used on reclaimed sites is an index devel¬ 
oped by Gale et al. (1991) for the assessment of habitat condi¬ 
tions for the cultivation of white spruce (Picea glciuca (Moench) 
Voss.) and an index described in the work of Burger et al. (1994) 
developed for the assessment of habitat conditions on reclaimed 
sites in the Appalachian Mountains for the cultivation of pine 
(Pinus virginiana Mill). 

In this work, the Mine Soil Quality Index (MSQI) was pro¬ 
posed for the assessment of initial mining soils and classification 
of habitats on post-mining sites afforested with a monoculture of 
the Scots pine (Pinus sylvestris L.). The MSQI may be useful as 
a tool for planning species composition, both for afforestation 
and pine monoculture conversion in the next generation of tree 
stands. 

Material and methods 

Study sites 

Field research was conducted on monoculture stands of the Scots 
pine (Pinus sylvestris L.), ranging from 12 to 30 years of age on 
the following reclaimed and afforested post-mining sites on dif¬ 
ferent parent rocks (substrate), including dominant types in Cen¬ 
tral Europe (Quaternary, Neogene and Carboniferous strata). 
First, a total of 32 square research plots measuring 10 m x 10 m 
each were set up and from there, four replications x 8 RMS sub¬ 
strate (parent rock) variants were established on four mine sites: 
(1) the hilltop of an external waste heap at Belchatow Lignite 
Mine; (2) a spoil heap of Smolnica hard coal mine; (3) the bot¬ 
tom of Szczakowa sand pit; and, 4) a waste heap of Piaseczno 
open-cast sulphur mine. Detailed site characteristics have been 
provided in Table 1. 

Soil study 

As part of the soil studies soil profiles on all study plots (total 32 
profiles) were exposed in pits to a depth of 110 cm. Additionally 
for composed soil sample preparation an extra 5 bore holes were 
made on each plot with a soil auger (Eijkelkamp set) on a grid. 
These composed soil samples were taken (1.0 kg mass of sample) 
to detennine basic soil properties at depths of 0-8 cm (Ai - ini¬ 
tial organic-mineral horizons); 8-50 cm and 50-110 cm (C - 
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parent materials/rocky spoils). Independently, soil bulk density 
(BD) was sampled by the core method, using standard sharpened 
steel cylinders (250 cm 3 ) (according to procedures in Dc Vos et 
al. 2005) with 3 replications for each layer (0-8; 8-50 and 
50-110 cm deep) in the RMS pits (one soil pit for each soil- 


substrate variant). In the lab, soil samples were dried and brushed 
through a 2.0 mm sieve. Samples of organic horizons OLf (anno¬ 
tated with FAO syntax where L = leaf; f = fermenta¬ 
tion/fragmentation humus layers) were collected after litterfall 
from 1 m xl m quadrates with 3 replications for each 100 nr. 


Table 1. Site characteristics 


Feature 

Belchatow spoil heap 

(BEL) 

Smolnica spoil heap 

(SMOL) 

Szczakowa sand pit 

(SCZ) 

Piaseczno spoil heap 

(PIAS) 

Location (GPS) 

N51 13.196 E19 25.569 

N50 15.078 E18 31.021 

N50 14.880 E19 23.932 

N50 33.622 E21 34.185 

Elevation (m.a.s.l.) 

350 

226 

282 

197 

Mean annual precipi- 

tation (mm) 

580 

702 

700 

650 

Mean annual tern- 

perature (°C) 

7.6 

7.7 

8.0 

7.0 


Parent material (sub¬ 
strate) 


BEL-QL: Quaternary loamy sands, 
loam, bouldery clay; 

BEL-NS: Neogene sands with loam, 
carbonated and sulphurised 


SMOL-CF and SMOL-NFC: 

Carboniferous deposits: mud¬ 
stones, sandstones and carbona¬ 
ceous shales 


SCZ-QS: fluvioglacial Quater¬ 
nary sands; 

SCZ-QLS: fluvioglacial Qua¬ 
ternary loamy sands 


PIAS-QS+NC: a mixture 
of Quaternary sands and 
Neogene Krakowiec Beds 
formation clays; 
PIAS-QLS: Quaternary 

loosjjsands_ 

PIAS-QS+NC and PIAS- 
QLS 

fertilisation (80 kg N, 50 kg 
P, 60 kg K ha 4 ); 2-years 
white legume cultivation 
(sowing 200 kg ha 4 ); plant¬ 
ing of pine trees. 


BEL-QL 


BEL-NS 


SMOL-CF 


SMOL-NFC SCZ-QLS and SCZ-QS 


Reclamation and 
reforestation treat¬ 
ments 


fertilization (60 kg fertilization and fertilization (100 unfertilized 

N, 70 kg P, 60 kg K afforestation as in kg N, 20 kg P, 20 substrates, 

ha 4 ); sowing the case of BEL-QL kg K ha 4 ); self- selfsown seed- 

grasses and legu- + neutralisation with sown seedlings of lings of pine, 

minous (60 kg ha* bog lime (60 Mg ha* pine. 

): planting of pine 1 ). 
trees. 


organic amendment addition 
(300 m 3 ha 4 organic horizons 
from forest soil, approx. 1.0 % 
C 0 rg content); liming (1.5 Mg 
dolomite ha 4 ); fertilization (140 
kg N, 130 kg P, 150 kg K ha* 1 ); 
2-years legume cultivation (240 

Ji 2 ilil 2 ; _P' il'll i 1 £ [222 ii+_ 


Pine tree stand age 
(years)_ 


17 


12 


30 


30 


21 


23 


30 


Next, mixed samples of organic horizon for each plot were 
prepared for laboratory testing. The basic soil parameters were 
determined in the samples using laboratory procedures: 

• particle size distribution was determined by hydrometer 
analysis and sand fractions by sieving; 

• the pH of the soil was measured with a combination 
electrode in suspensions of 1.0 mol L* 1 KC1 (pH K ci) (1:2.5 
mass/volume ratio) after 24 hour equilibration; 

• organic soil carbon (C org ) and total nitrogen (N t ) were 
assayed with a Leco CNS 2000 Elemental Analyzer; and 

• basic exchangeable cations (Na + , K + , Ca 2 ~, Mg 2+ ) were 
extracted with 1 mol L _1 NH 4 OAc. 

Samples were then mixed with a small portion of extractant 
and equilibrated. After 24 hours, the suspensions were filtered, 
the soils were washed with additional extractant, and the total 
volume was made up to 100 mL (Jackson, 1958). The concentra¬ 
tion of cations was determined by atomic absorption spectros¬ 
copy (AAS) with a Varian Spectrophotometer. Phosphorus in the 
fonn available to plants (P av ) was assayed using the Egner- 
Riehm method in calcium lactate extract ((CHgCHOHCOOl^Ca) 
acidified with hydrochloric acid to pH 3.6 and using the colori¬ 


metric method (Ostrowska et al. 1991) with VARIAN Inc., Cary 
300 UV-Vis Spectrophotometers. 

Plant community studies 

Vegetation coverage on the study plots was determined using the 
Braun-Blanquette method (100 m 2 phytosociological surveys on 
each experimental plot, total of 32) (Wikum and Shanholtzer 
1978). Next, ecological indicators of moisture (M) and fertility 
(F) for the forest floor species composition were calculated based 
on ecological numbers of vascular plants in Europe (Ellenberg 
2009). In the analyzed cases, an 5-grade scale was used where 
the intensity of particular factors increases from (0) 1 to 5 (6) 
(Zarzycki et al. 2002). 

As is widely known many plant species are good indicators of 
habitat conditions. Detailed studies conducted 60 years ago 
allowed tocalibrate many plant species with reference to habitat 
conditions; Ecological indicator values for Central Europe 
(Ellenberg 2009) have been used successfully by many authors 
(Dzwonko 2001). In most cases, a 5-grade scale was used where 
the intensity of particular factors increases from (0) 1 to 5 (6) 
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(Zarzycki et al. 2002). Here, to calculate the values of the given 
indices, the species abundance index according the Braun- 
Blanquette scale was converted into numerical values: cover 
above 75% corresponds to "5", cover range from 50%-75% to 
"4", 25%-50% to "3"; very abundant or 5%-25% to "2"; 
numerous individuals but not abundant to "1"; rare and not 
abundant to "0.5"; few individuals to "0.2" and very rare (one) 
corresponds to "0.1". Calculations were made as a weighted 
average of ecological indicator values, first for an individual 
phytosociological survey on each plot and later for mean values 
for substrate type (variants) on the sites. 

Forest floor vegetation biomass was determined using the 
harvest method on 96 sample sub-plots (1 m 2 x 3 replications 
distributed along the diagonal of each of 32 soil study plots 
measured 10 m x 10 m) in the middle of the vegetation season 
(July). Mixed samples of vegetation (aboveground biomass of 
herbaceous and shrubs, Biom H+S) were collected to determine 
the water content (%) and calculate dry mass in the laboratory. 

Statistical procedures 

Data were analyzed with Statistica 9.1 software (StatSoft Inc. 
Software, 2009). Significant differences between the mean val¬ 
ues of basic soil characteristics and MSQI from differing groups 
of soil-substrate (parent rock) variants were tested by an 
ANOVA test, preceded by a Shapiro-Wilk test of normality, and 
Levene’s test of variance homogeneity. The ANOVA test was 
followed by multiple pairwise comparisons using Tukey's HSD 
(honestly significant difference) post-hoc test. Based on the re¬ 
sults of Tukey's HSD test homogeneous subsets in tenns of index 
values MSQI were distinguished. Correlation analysis between 
MSQI and plant features for the validation model was performed 
using the Pearson’s correlation coefficient r (at the significance 
level p=0.05). 

The development and calculation of mine soil quality index 
(MSQI) 

Following the studies and concepts presented by Gale et al. 
(1991), Burger and Kelting (1999), Brozek et al. (2011) some 
assumptions were made to develop the MSQI index. The com¬ 
ponents of the soil assessment index include basic soil properties 
contributing to its fertility such as soil texture, nutrient availabil¬ 
ity, acidity converted into volumetric units, and sub-indices. 
Furthermore, each of the sub-indices was weighted with consid¬ 
eration of its estimated impact on the final assessment and stan¬ 
dardised according to absolute values ranging from 0.0 to 1.0. 
Each feature was standardised using a linear function, taking into 
account three approaches: first, more (increase in feature value) 
is worse; second, more of a given index is better for the plant; 
third, at first more is better, then there is an optimum, and further 
increase is worse. 

MSQI takes account of the following features (Table 2): the 
stock of soil texture fractions; nutrient availability, acidity (pH) 
and biological activity (expressed by C org -to-N t ratio). The value 
of each feature was first defined for particular horizons, i.e.: 0-8 
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cm (organic mineral initial horizons Ai); 8-50 cm (parent rock C 
horizon, a layer most intensively penetrated by tree root systems 
in mine soils); and 50-110 cm (parent rock C horizon, approxi¬ 
mate root range). The index referring indirectly to the decay rate 
of precipitate in the form of organic matter and soil biological 
activity (i.e., the C org -to-N t ratio), was defined only for the forest 
litter and partially decomposed humus layers (OLf). 

The developed MSQI index is the sum of sub-indices 
weighted with a relevant validity coefficient factor (wt) provided 
for a given feature, where the maximum value of 1.0 corresponds 
to 100%. 

Finally, a general equation was applied to calculate the MSQI 
(Equation 1): 

MSQI= (1st wt(0.3)) + (Ina wt(03)) + (IP wt(0.1)) + 

(lac wt(0.2)) + (Iba ivt(O.l)) (1) 

where: 

wt is weight as a validity coefficient factor for individual sub-index; 
1st is sub-index of soil texture, which is the sum of Isi, the silt sized 
fraction stock (SFS, 0.05-0.002-mm fraction) and Icl, the clay sized 
fraction stock (CFr, <0.002-mm fraction) calculated according to the 
total stock of these fractions in soil up to 110 cm profile deep 
(Mg-ha~'T 10-cm 4 ). However, the value of this sub-index is reduced 
by the Is skeletal sub-index defined according to the % of coarse 
fragment (CFr, >2.0 mm) components (Equation 2): 

1st = (Isi + Icl) - Is (2) 

however, 1st >0 because with extreme values (small share of silt and 
clay sized fraction with very high coarse fragment CFr share, i.e., > 
35 % ) the 1st sub-index would be negative and so in such cases 1st is 
assumed to be 0.0, wt = 0.3; 

Ina - nutrient availability index refers to the nutrient content in soil 
volume (ECS, exchangeable cation form stock Ca 2 ', Mg 2+ , K , Na + 
up to 110 cm soil profile deep), wt = 0.3; IP available phosphorus 
sub-index based on Available Phosphorus Stock (APS) in soil vol¬ 
ume (Mg-ha ' l 10-cm' 1 ), wt = 0.1; lac - acidity sub-index based on 
hydrogen stock (H + S) in soil volume (in units in kg-ha ''II 0-cm 1 ); 
however H + S has been calculated based on pH K ci, which constitutes a 
measure of hydrogen ions in soil water solution ratio 1:2.5. Such a 
manner of determining acidity allowed the researchers to avoid tak¬ 
ing the arithmetic mean of pH from soil horizons (pH is a logarithmic 
value); wt = 0.2; Iba - soil biological activity index based on C 0rg -to- 
N, ratio in OLf horizon, wt 0.1: 

To automate the calculation of MSQI based on a scale of 
valuation, the IF function was used to perfonn logic tests on 
values and formulas in an Excel MS Office (2007) spreadsheet. 

Evaluation of individual MSQI variants at different wt for sub¬ 
indices was validated by the best correlations with plant commu¬ 
nity characteristics affected by natural succession, including 
aboveground biomass of forest floor plant and ecological indica¬ 
tor of vascular plants (moisture M and fertility F). In this work, 
the tree stand characteristics were not considered for validation 
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of MSQI, because they heavily depended on the age of trees, 
which varied from 12 to 30 years at the post-mining site (Table 
1 ). 

The best results (at a significance level of p= 0.05) were at 
“wF value respectively: 0.3 for 1st and Ina, 0.2 for lac, and 0.1 
for IP and Iba. In this case, the values of MSQI correlated sig¬ 


nificantly with the ecological indicators: moisture (M) (r =0.39) 
and fertility (F) (r =0.42) (Fig. 1 a, b). The values of MSQI cor¬ 
related significantly (at p =0.05), as well, with aboveground 
biomass of forest floor species (herbaceous and shrub vegetation, 
Biom H+S, r=0.39) (Fig. 1, c). 



Fig. 1 . Correlations between mine soil quality index MSQI and ecological indicators of Fertility F (a), Moisture M (b) and aboveground forest floor spe¬ 
cies community biomass (herbaceous and shrub vegetation Biom H+S) (c) in a validation model. According to Pearson’s test with n = 32, at p =0.05 the 
statistical correlation at 0.3 = < r xy< 0.5 is average. 


Table 2. Scale and list of individual partial index ranges which constitute the mine soil quality evaluation according to MSQI formula. 


Group of soil properties and sub-indices of MSQI 


Texture - stock of particle size fractions: 

: 1st = 

(Isi + Icl) - Is 



Nutrient stock 


Acidity 


Biological activity 

SFS 1 

(0.05-0.002 mm) 

Mg-ha'M 10-cm' 1 

Isi 

C1FS 

(<0.002 mm) 
Mg-ha'M 10-cm' 1 

Icl 

CFr 

(>2.0 mm) 

% 

Is 

ECS 

Mgha’ 1 -! 10-cnT 

j lab 

APS 

Mg-ha" 1 T 10-cm’ 1 

IP 

H + S 

kg-ha _1 T 10-cm’ 1 

lac 

Corg - tO - Nt 

Iba 

<250 

0.05 

<300 

0.05 

<15.0 

0 

<5.0 

0.15 

<0.005 

0.15 

>20.0 

0.2 

<28.0 

i 

251-600 

0.1 

301-600 

0.1 

15.0-25.0 

0.1 

5.1-10.0 

0.3 

0.0051-0.0100 

0.3 

10.1-20.0 

0.3 

28.1-35.0 

0.9 

601-900 

0.15 

601-1000 

0.15 

26.0-35.0 

0.2 

10.1-15.0 

0.45 

0.0101-0.0200 

0.45 

2.1-10.0 

0.4 

35.1^10.0 

0.6 

901-1400 

0.2 

1001-1500 

0.2 

>35.0 

0.3 

15.1-20.0 

0.6 

0.0201-0.0360 

0.6 

1.1-2.0 

0.5 

40.1-50.0 

0.5 

1401-1600 

0.25 

1501-2000 

0.25 



20.1-35.0 

0.75 

0.0361-0.0600 

0.75 

0.5-1.0 

1 

50.1-60.0 

0.4 

1601-2500 

0.3 

2001-3000 

0.3 



35.1-50 

0.9 

0.0601-0.1000 

0.9 

0.002-0.5 

0.6 

60.1-70.0 

0.2 

2501-3700 

0.35 

3001-3500 

0.4 



>50 

1 

>0.1000 

1 

<0.002 

0.3 

>70.0 

0.1 

3701^1300 

0.4 

3501—4000 

0.5 











4301-8500 

0.45 

4001-8500 

0.2 











>8500 

0.5 

>8500 

0.1 












1 SFS - Silt Sized (0.05-0.002 mm) Fraction Stock; C1FS — Clay Sized (<0.002 mm) Fraction Stock; CFr — Coarse Fragments (>2.0mm); ECS — Exchangeable 
Cations Form Stock; APS — Available Phosphorus Stock; H + S — Hydrogen Stock; 


Results 

Plant characteristics 

The mean plant ground cover abundance was varied and ranged 
from 2 (SCZ-QS) to 56 % (SMOL-CF); the mean number of 
undergrowth species of vascular plants ranged from 13 (PIAS- 
QLS) to 45 (SMOL-NCF). In total, 131 plant species were iden¬ 
tified in the post-mining sites, 29 of which were classified as 


forest species, 80 as ruderal and 22 as grassland, respectively. 
The dominant species in the communities were: Calamagrostis 
epigejos, Cardaminopsis arenosa, Cirsium arvense, Dactylis 
glomerata, Deschampsia flexuosa, Festuca ovina, Festuca rubra, 
Hieracium pilosella, Poa compressa, Tussilago farfara, Vaccin- 
ium myrtillus respectively depending on the substrate (Table 3). 
However, the remaining species had no clear phytosociological 
relation or were characteristic for various other vegetation types. 
The value of the ecological indicator of moisture (M) ranged 
from 2.60 (BEL-NS) to 3.18 (PIAS-QS+NC) and the value of the 
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fertility indicator (F) ranged from 2.35 (SMOL-CF) to 3.43 biomass (Biom H + S) ranged from 0.13 (SCZ-QS) to 0.345 
(PIAS-QS+NC) (Table 3). The forest floor plant aboveground (Mg-ha' 2 ) (SMOL-CF) (Table 3). 


Table 3. Selected plant community features 


Feature of plant community 

Plant ground 

cover abun- 

Number of 

vascular plant 

Dominant species 1 

Aboveground 

biomass (forest 

Moisture 

(M) 3 

Fertility 

(F) 



dance %) 

species (pcs.) 


floor) (Mg-ha' 1 ) 

The BELCHATOW 

BEL-QL 

8 

37 

Dactylis glomerata; Tussilago farfara; 

Cirsium arvense; Deschampsia flexuosa 

0.135 (0.091) 2 

2.98 (0.10) 4 

3.35 (0.06) 

spoil heap 

BEL-NS 

5 

17 

Festuca rubra; Hieracium pilosella; Poa 

compressa 

0.033 (0.021) 

2.60 (0.20) 

2.83 (0.34) 

The SMOLNICA 

SMOL-CF 

56 

30 

Deschampsia flexuosa; Festuca ovina; 

Vaccinium myrtillus 

0.345 (0.241) 

3.03 (0.05) 

2.35 (0.13) 

spoil heap 

SMOL-NCF 

26 

45 

Deschampsia flexuosa; Calamagrostis 
epigejos; Festuca ovina 

0.153 (0.102) 

3.05 (0.13) 

2.53 (0.21) 

The SZCZAKOWA 

sand pit 

SCZ-QLS 

SCZ-QS 

4 

2 

34 

22 

Deschampsia flexuosa; Hieracium pilosella 

Cardaminopsis arenosa; Deschampsia 

flexuosa; Hieracium pilosella 

0.015 (0.010) 

0.013 (0.011) 

2.90 (0.22) 

2.63 (0.19) 

2.98 (0.17) 

2.58 (0.29) 

The PIASECZNO 

PIAS-QS+NC 

34 

34 

Solidago gigantea; Festuca rubra; 

Equisetum arvense 

0.170 (0.245) 

3.18 (0.10) 

3.43 

(0.10) a 

spoil heap 

PIAS-QLS 

4 

13 

Festuca rubra 

0.073 (0.062) 

3.10(0.12) 

3.18 

(0.05) b 


1 dominat species in phytosociological surveys on plots represent site variant; 


2 0.135 (0.091) - mean and standard deviation of aboveground forest floor biomass (dry mass) (Mg-ha 4 ) 

2 ecological indicators of moisture (M) and fertility (F) based on Ellenberg (2009) (description in methods section); 

4 2.98 (0.10) - mean and standard deviation (number of replications for phytosociological survey on plots represent site variant n = 4) 


Basic soil properties 

The properties of the investigated RMS are rather varied 
(statistical differences of basic soil parameters are given in Table 
4). The percentage of silt-sized fraction (0.05-0.002 mm) ranged 
from 2.0% to 36.0% and clay sized was (<0.002 mm) between 
1.0 and 25.0%. Organic carbon (C org ) and nitrogen N t contents in 
the Ai horizon (0-8 cm depth) varied from 2.3 to 166.1 g-kg' 1 and 
0.15 to 4.04 g-kg' 1 , respectively and C org -to-N t ratio in the OLf 
horizon between 29.3 and 83.6. Soil pH in 1 M KC1 ranged from 
3.13 to 7.50. The samples also differed considerably in terns of 
their cation exchange properties. The total exchangeable bases 
(TEB) were from 1.20 to 27.41 cmol^ykg' 1 , cation exchangeable 
capacity (CEC) from 2.18 to 28.0 cmol^-kg" 1 and base saturation 
(BS) from 25.96 to 98.14%. The content of available phosphorus 
(P av ) ranged from 0.10 to 2.68 g-kg' 1 (Table 4). 

Value of transfonned soil properties for MSQI calculation 

Depending on the silt fraction content in the horizons of the in¬ 
vestigated soils, the average stock of this fraction (SFS) ranged 
from 343 (SCZ-QS) to over 4827 Mg-ha'-l 10-cm' 1 (BEL-QL) 
(Table 5). The clay fraction stock (C1FS) in the investigated soils 
ranged from 263 (SCZ-QS) to 2775 Mg-ha' 1 -110-cm' 1 (SMOL- 
CF) of the profile (Table 5). For the investigated RMS, the mean 
value of Isi ranged from 0.10 to 0.45 and Icl stock ranged from 
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0.05 to 0.5 (Table 5). The coarse fragment fraction CFr (frag¬ 
ments of rocks, stones, gravel of >2.0 mm in diameter) in RMS 
ranged from 0.0% (in most soils) to 80.0% (SMOL-CF). Eventu¬ 
ally the mean value of 1st index (according to equations 2) 
ranged from 0.20 (SCZ-QS and PIAS-QLS) to 0.60 (BEL-QL) 
(Table 5). Exchangeable cation fonn stock (ECS) include: cal¬ 
cium Ca 2+ , magnesium Mg 2+ , potassium K + and sodium Na~ of 
investigated soil ranged from 2.8 (SZC-QS) to 99.2 Mg-ha' 1 -110- 
cm' 1 (BEL-QL) (Table 5). 

In the investigated RMS, as in most natural soils, calcium 
(Ca 2+ ) was the predominant nutrient. Available phosphorus stock 
(APS) in the investigated mining soils ranged from 0.0025 
Mg-ha' 1 -! 10-cm' 1 in soils forming on Carboniferous unfertilized 
sediments (SMOL-NCF) to 0.4800 Mg-ha' 1 -110-cm' 1 in soils 
forming on acidic, sulphurous neogene sands, neutralised by bog 
lime (BEL-NS) (Table 4). The mean value of IP in investigated 
RMS ranged from 0.30 (for soils in SMOL-NCF) to 1.0 (BEL- 
QL and BEL-NS) (Table 5). 

The acidity of the investigated soils expressed by hydrogen 
stock in terms of volume (H + S) ranged on average from 0.0013 
(BE-QL) to over 43.5250 kg-ha‘’-l 10-cm' 1 (Table 5). The mean 
lac for investigated soil substrate ranged from 0.2 (SMOL-CF 
and SMOL-NCF) to 1.0 (SCZ-QLS and SCZ-QS). The mean 
value of C 0rg -to-N, ratio in OLf horizons of the investigated soils 
ranged from 31.5 (SMOL-CF and SMOL-NCF) to 83.6 (BEL- 
NS) (Table 2), and Iba sub-index for investigated soils ranged 
respectively from 0.1 (BEL-QL and BEL-NS) to 0.9 (SMOL-CF; 
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SMOL-NCF; PIAS-QS+NC and PIAS-QLS) (Table 5). 


Table 4. Basic soil characteristics 


Feature 

BD 

(gem 

f 3 ) 

Corse frag¬ 
ments (%) 

Silt sized fraction 

(0.05-0.002 mm) (%) 

Clay sized fraction 

(<0.002 mm) (%) 

pHKCl 

C org (g kg 'j 


horizon 

Ai 

c 

Ai, C 

Ai 

C 

Ai 

C 

Ai 

c 

OLf 

Ai 

BEL-QL 

1.64 a ± 0.06 

1.69 a ± 0.07 

0 

32 b ± 11 

30 a ± 8 

3 a ± 3 

5 a ±4 

7.36 c ± 0.19 

7.50 c ± 0.12 

478.5 C ± 5.5 

6.1 a 

±2.8 

BEL-NS 

1.62 a ± 0.06 

1.71 a ± 0.07 

0 

8 a ± 2 

17 a ± 19 

3 a ± 2 

6 a ± 2 

4.94 ac ± 1.75 

4.69 b ± 1.50 

407.5 a ± 53.4 

2.T 

±0.4 

SMOL-CF 

1.00 a ± 0.06 

1.1 l a ± 0.02 

70±80 

36 b ±2 

33 a ± 2 

24 c ± 2 

25 a ± 2 

3.40 ab ± 0.20 

3.60 a ± 0.43 

436.0 a ± 35.0 

166.1 

a ± 13.0 

SMOL-NCF 

1.07 a ± 0.12 

1.04 a ± 0.04 

70±80 

36 b ±2 

31 a ±4 

24 c ± 2 

25 a ± 2 

3.40 ab ± 0.20 

3.13 a ± 0.21 

436.0 a ± 35.0 

166.1 

a ± 13.0 

SCZ-QLS 

1.54 a ± 0.13 

1.58 a ± 0.02 

0±5 

8 a ± 3 

12 a ± 8 

4 a ± 1 

5 a ± 3 

4.15 ab ± 0.06 

4.60 ab ± 0.17 

316.6 a ± 96.2 

4.6 a 

± 1.3 

SCZ-QS 

1.63 a ± 0.06 

1.52 a ± 0.11 

0±5 

5 a ± 1 

2 a ± 1 

l a ± 1 

2 a ± 1 

5.91 ac ± 2.05 

e.lS^i 1.61 

333.3 a ± 72.3 

2.3 a 

± 1.2 

PIAS-QS+NC 

1.41 a ± 0.16 

1.29 a ± 0.10 

0+10 

1 l a ± 6 

8 + 4 

9 b ± 2 

7 a ± 2 

5.56 ac ± 1.36 

6.40 hc ±0.78 

297.8 ab ± 40.8 

13.2 

a ± 1.9 

PIAS-QLS 

1.53 a ± 0.06 

1.62 a ± 0.01 

0+5 

4 a ± 2 

3+ 1 

4 a ± 1 

3 a ± 1 

5.39 ac ± 1.60 

e^iO^O 

329.0 a ± 103.4 

4.7 a 

±0.4 

Feature 

N, (g-kg" 1 ) 

C ore -to-N t 

TEB (cmol c kg ') 

CEC (cmol c kg'') 

BS% 

P„(gkg l ) 


horizon 

OLf 

Ai 

OLf 

Ai 

C 

Ai 

C 

Ai 

c 

Ai 


C 

BEL-QL 

5.79 ab ± 0.32 

0.53 a ± 0.03 

82.8 a ± 3.8 

27.41 b ± 5.68 

27.18 c ± 4.16 

28.00 c ± 5.68 

27.68 c ± 4.23 

97.86 b ± 0.42 

98.14 d ±0.14 

2.68 b ± 1.62 

1.55 a 

± 1.20 

BEL-NS 

4.87 a ± 0.48 

0.15 a ± 0.02 

83.6 a ± 5.7 

4.07 a ± 2.00 

3.28 ab ± 0.97 

5.69 a ± 1.39 

5.49 a ± 0.85 

68.81 ab ± 20.31 

61.02 bo ±22.72 

1.93 b ± 0.94 

2.66 a 

±0.85 

SMOL-CF 

13.84 c ± 0.42 

4.04 a ± 0.38 

31.5 a ± 2.6 

5.15 a ± 1.82 

10.62 b ± 2.74 

19.40 bc ± 6.49 

21.30 bc ± 2.95 

28.37 a ± 13.12 

47.05 ab ± 10.39 

0.18 a ± 0.16 

0.21 a 

±0.22 

SMOL-NCF 

13.84 c ± 0.42 

4.04 a ± 0.38 

31.5 a ± 2.6 

5.15 a ± 1.82 

4.53 ab ± 1.67 

19.40 bc ± 6.49 

16.24 b ± 4.69 

28.37 a ± 13.12 

25.96 a ±4.91 

0.18 a ± 0.16 

0.21 a 

±0.22 

SCZ-QLS 

7.70 b ± 2.15 

0.40 a ± 0.06 

41.0 a ± 5.2 

1.20 a ± 0.37 

2.10 a ± 0.32 

2.21 a ± 0.72 

3.13 a ± 0.49 

56.68 ab ± 18.81 

66.72 bc ±8.06 

0.12 a ± 0.03 

0.1 l a 

±0.03 

SCZ-QS 

7.83 b ± 1.75 

0.30 a ± 0.09 

42.8 a ± 5.1 

2.59 a ± 2.72 

1.28 a ± 0.92 

3.68 a ± 2.44 

2.18 a ± 1.11 

52.05 ab ± 38.49 

49.63 ab ± 12.12 

0.24 a ± 0.21 

0.1 o a 

±0.03 

PIAS-QS+NC 

10.14 bc ± 1.11 

0.86 a ± 0.11 

29.3 a ± 1.2 

11.29 a ± 9.60 

10.2 l b ± 7.95 

14.10 ab ± 8.36 

11.75 ba ± 7.47 

73.49 ab ± 16.51 

83.34“* ±9.31 

0.82 a ± 0.83 

0.43 a 

±0.23 

PIAS-QLS 

9.44 bc ± 2.96 

0.35 a ± 0.05 

35.0 a ± 2.9 

3.79 a ± 3.77 

2.69 ab ± 1.53 

5.44 a ± 3.47 

3.73 a ± 1.36 

59.67 ab ± 21.95 

68.01 c ± 12.97 

0.17 a ± 0.23 

0.16 a 

±0.08 


Notes: OLf- initial organic horizon; Ai - initial organic - mineral horizon at 0-8 cm deep; C - parent rock horizon, mean values for 8—50 and 50-110 cm layers 
cm deep; 1.64' ± 0.06 - sample mean and standard deviation ( n 4, for BD n=3 for each horizon), for coarse fragment range of % content in volume of soil; site and 
substrate variants abbreviations — see description in the Table 1; Means of soil-substrate features in horizons with the same letter are not significantly different at p= 
0.05 by Tukey test. 


Table 5. Values of the sub-indices of mine soil quality index (MSQI) calculated on the basis of mine soil characteristics. 


Mine sites 

and substrate 

CFr 

% 

SFS 1 

C1FS 

Is 

Isi 

Icl 

ECS 

APS 

Ina 

IP 

Corg - tO - Nt 

Iba 

H + S 

lac 

Mg-ha' 1 • 110-cm" 1 

kg-haM 10-cm' 1 

Mg-ha' 1 

•110-cm* 1 

BEL-QL 

0 

4827 2 ± 1094 

771 ±87 

0.00 

0.45 

0.15 

99.2 ± 13.4 

0.2890 ±0.2128 

1.00 

1.00 

82.8 ±3.8 

0.1 

0.0013 ±0.0005 

0.30 

BEL-NS 

0 

2418 ± 1145 

982 ± 397 

0.00 

0.40 

0.20 

10.3 ±3.0 

0.4800 ±0.1600 

0.45 

1.00 

83.6 ±5.7 

0.1 

27.1510 ± 40.8169 

0.50 

SMOL-CF 

80 

3994 ±111 

2775 ±312 

0.30 

0.40 

0.40 

29.3 ±7.6 

0.0235 ± 0.0230 

0.75 

0.45 

31.5 ±2.6 

0.9 

13.8250 ± 13.4441 

0.20 

SMOL-NCF 

80 

3717 ± 239 

2563±172 

0.30 

0.40 

0.30 

9.1 ±4.0 

0.0025 ± 0.0038 

0.45 

0.30 

31.5 ±2.6 

0.9 

43.5250 ± 17.9632 

0.20 

SCZ-QLS 

0 

1952 ±700 

891±183 

0.00 

0.35 

0.15 

8.2 ± 1.2 

0.0182 ±0.0040 

0.30 

0.60 

41.0 ± 5.1 

0.5 

1.0450 ±0.6087 

1.00 

SCZ-QS 

0 

343 ±115 

263 ± 44 

0.00 

0.10 

0.10 

2.8 ± 1.0 

0.0213 ±0.0074 

0.15 

0.45 

42.9 ±5.1 

0.5 

0.3388 ±0.3909 

1.00 

PIAS-QS+NC 

5 

1172 ±706 

1001 ±410 

0.00 

0.30 

0.20 

27.4 ±21.2 

0.0657 ±0.0367 

0.75 

0.45 

29.3 ± 1.3 

0.9 

0.0908 ±0.1510 

0.60 

PIAS-QLS 

2 

387±102 

556 ±90 

0.00 

0.10 

0.10 

8.3 ±2.4 

0.0302 ±0.0125 

0.30 

0.45 

35.0 ± 2.8 

0.9 

0.2183 ±0.2578 

0.60 


1 SFS - Silt Sized (0.05-0.002 mm) Fraction Stock; C1FS - Clay Sized (<0.002 mm) Fraction Stock; CFr - Coarse Fragments (>2mm); ECS - Exchangeable 
Cations Form Stock; APS - Available Phosphorus Stock; C or g-to-N t - organic Carbon to total Nitrogen ratio; H S - Flydrogen (H ) Stock; 

2 Values of transformed soil characteristics for site and substrate variants are mean with standard deviation (n= 4). 


Discussion 

As one can see from the above characteristics, comparison of soil 
quality using particle soil properties is difficult. It is an ex¬ 
tremely important to use soil quality indices to make objective 
comparisons of mine soils developed on varied parent rock mate¬ 


rial and at different post-mining facilities such as external spoil 
heaps (in sulphur, lignite and hard coal mining) and sand open 
castpits. Deposits that constitute parent rock in soils on post¬ 
mining sites are formed from Quaternary, Neogene and Carbon¬ 
iferous deposits with different physical and chemical features. 
When numerically assessing various features of transformed 
soils into sub-indices, what was taken into account was the scale 
of variability range of these features in the analysed data set. 
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What was also taken into account were the requirements and 
ecological properties of the tested species, the Scots pine (Pinus 
sylvestris L.), which dominates the afforested post-mining sites 
in this part of Europe. The species' wide ecological amplitude, its 
pioneering in ecological succession and good adaptation poten¬ 
tial for difficult habitat conditions (Farjon 2005) are the reasons 
for its widespread use in afforestation of post-mining sites (Rno- 
che 2005; Baumann et al. 2006; Pietrzykowski 2010; 
Pietrzykowski and Socha 2011). 

Soil texture is one of the most important criteria for the classi¬ 
fication of soil deposits. This feature is closely connected and 
decisive for other soil properties. Among the various fractions, 
the most important role is played by silt and clay. Silt fraction 
(0.05-0.002 mm) in sands increases water capacity and capillary 
conduction, whereas in clays it reduces swelling, viscosity and 
plasticity. Colloidal clay admixture increases cohesion and plas¬ 
ticity and reduces percolation and water-permeability. Excessive 
content of the clay fraction, however, adversely affects the air 
and water properties of soils and water availability for plants. 

When developing the soil quality evaluation scale, the silt in¬ 
dex (Isi) was based on the assumption that the larger stock of this 
fraction is “the better”. When developing the soil quality evalua¬ 
tion scale for the clay index (Icl), increases in the content and 
stock of the clay fraction was considered favourable only up to a 
certain limit. It was assumed that it was most beneficial from 
3501-4000 Mg-ha‘'-l 10-cm' 1 , and over this range the Icl value 
fell. If the clay content in profile exceeds 8500 Mg-ha' 1 -110-cm' 1 
cm (that is if the clay fraction constituted on average over 60 % 
in bulk density of the investigated soils), the Icl sub-index value 
will drop to 0.2. The coarse fragment fraction CFr (fragments of 
rocks, stones, gravel of >2.0 mm in diameter) in RMS decreased 
the value of soil texture index (1st). In determining the partial 
value of the Is, coarse fraction was not calculated in advance, as 
in the case of silt and clay fractions, but the valuation was made 
on the basis of the percentage (% in sample volume). If the con¬ 
tent of silt and clay fraction was very small and the CFr very 
high in tipped mine deposits, the 1st value in the valuation of soil 
texture could not be negative (1st >0), and in such cases it is as¬ 
sumed that 1st = 0.0. 

Another feature included in the valuation of mine soil was 
availability of nutrients for plants. This feature obviously has a 
direct impact on growth and nutrition of the introduced tree 
stands. Nutrient availability sub-index Ina refers to exchangeable 
cation form stock (ECS): calcium Ca 2+ , magnesium Mg 2+ , potas¬ 
sium K + and sodium Na + . The available phosphorus (P av ) is also 
included in the group of basic macronutrients which are fre¬ 
quently deficient in mine soils (Heinsdorf 1996; Daniels and 
Stewart 2000; Pietrzykowski 2010). In natural soils, phosphorus 
occurs in mineral form and its source is mainly fluoroapatite; in 
organic form, it is connected with plant and animal remains. 

In initial mine soils, little organic matter has accumulated and 
apart from its natural content in rock, which forms spoil heaps 
and excavations, some amounts of phosphorus are provided in 
the form of mineral fertilising. The impact of fertilising tree- 
stands with phosphate may last from several to about a dozen 
years (Baule and Fricker 1970, Heinsdorf 1996). The scale for IP 
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was estimated assuming that the increase of the index value is 
linear, that is “the more phosphorus, the better for the plants.” In 
some cases on sulphurous neogene sands (BEL-NS), neutralisa¬ 
tion by bog lime improved the IP sub-index values. Bog lime 
used for neutralization is a rock containing phosphorus bound 
with organic residue formed in the course of sedimentation in 
neogene lakes. 

Soil acidity expressed by pH is a feature very closely con¬ 
nected with other soil properties. The pH in 1MKC1 determina¬ 
tion is based on an assumption that hydrogen and aluminium 
ions go into solution in reaction with neutral salt (KC1), and the 
obtained value is quite stable in the vegetation season and less 
dependent on soil humidity, for example. This is of particular 
importance in the case of mine soils that frequently exhibit sig¬ 
nificant pH variability in vertical profiles. Forest tree species 
have different requirements regarding pH for both the optimum 
range and the tolerated range (ecological amplitude). 

The investigated species, the Scots pine, has a wide range of 
tolerance for pH variation, even within the soil profile (Farjon 
2005). In the valuation of an lac acidity index, it was assumed 
that the smaller the accumulation, the better, but only to a certain 
limit and then the lac index grew, and subsequently, below 1.0 
kg H T ha" 1 110 cm' 1 , the lac index dropped (Table 5). This as¬ 
sumption stems from the fact that most forest species grow in 
natural conditions in soils that are at least slightly acidic and, in 
the case of mine soils, both extremely low acidity (e.g. mining 
drainage) and alkalinity are adverse phenomena. 

An important criterion in the evaluation of RMS development 
in post-mining areas is the extent and dynamics of initial organic 
horizon fonnation. One measure of development of these hori¬ 
zons is the accumulation of soil organic Carbon (C org ) and total 
nitrogen (N t ) and the C org -to-N t ratio (Anderson 1977; Wali and 
Freeman 1973; Prosser and Roseby 1995; Li and Daniels 1994, 
Rumpel et al. 1999, Pietrzykowski 2008). The fresh litter horizon 
(OL) is clearly developed and the raw humus horizon (Of) is at 
an early development stage and hence the symbol OLf (fresh 
litter and raw humus horizon). 

The C or g-to-N t ratio may be interpreted as an indirect index of 
changes in the developed soil environment, including the inten¬ 
sity of soil organic matter transformation processes (Janssen 
1996). With high C„ rg -to-N, ratio, the stock of mineralised nitro¬ 
gen may be used by soil micro-organisms, and consequently the 
accumulation of nitrogen may be stopped. In the sediments 
freshly deposited on spoil heaps, the phenomenon of nitrogen 
deficiency (too high C org -to-N t ratio) may occur periodically 
(Wali 1999). In the case of RMS, wide C org -to-N t ratio is con¬ 
nected with slower decomposition of biomass produced by pio¬ 
neering plant communities (Schafer and Nielsen 1979). 

In the case of initial mine soils, there is however a difficulty in 
interpreting C org -to-N t ratio due to the occurrence of origin 
(“geogenic”) carbon that significantly alters the total balance of 
carbon and soil nitrogen (Chabbi et al., 2008). And so the pro¬ 
posed MSQI includes the C on ,-to-N t ratio in OLf horizon in the 
Iba partial index, which refers to “biological activity”. From the 
literature (Farjon 2005) it follows that in habitats suitable for the 
investigated species (Scots pine), the C org -to-N t ratio in litter is 
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on average above 60. Series of data indicate, however, that in 
optimum habitats and with a good supply of nitrogen, the C org -to- 
N, ratio in the Scots pine may be lower than 30 and vary dynami¬ 
cally with the decomposition of litter. Usually the amount of 
nitrogen in comparison to organic carbon increases in the course 
of litter decomposition. 

MSQI relations with features of vegetation 

As mentioned above, validation of soil quality indices and suit¬ 
ability in habitat classification may be done through correlation 
with plant community features (Warkentin 1995; Burger and 
Kelting 1999; Schoenholtz et al. 2000). It is well known that 
many plant species from succession are good indicators of habi¬ 
tat conditions. Detailed studies conducted 50 years ago allowed 
to calibrate many plant species with reference to habitat condi¬ 
tions and ecological indicator values for Central Europe (Ellen- 
berg first in 1965 and reprinted in Ellenberg 2009) were used 
successfully. The values of MSQI correlated significantly with 
ecological indicator values of vascular plants: moisture (M) and 
fertility (F). Variants on Carboniferous shales (SMOL-CF and 
SMOL-NCF), where the values of MSQI were higher while soil 
fertility (F) values were the lowest, had significant impact on the 
reduction of the correlation coefficient values between MSQI 
and ecological indicator of F. This indicates that predictions of 
site conditions based on plant community features in some cases 
may be lower than the assessment of potential soil quality. In 
such cases, it will be necessary to better plan potential habitat 
conditions for the introduced tree stands. 

The biomass of plant communities and the rate of production 
of biomass per unit area are other factors important for soil qual¬ 
ity and productivity validation. The quantity and spatial distribu¬ 
tion of biomass layers and plant constituents is of crucial signifi¬ 
cance in describing ecosystem productivity (Krebs, 1994). How¬ 
ever, the plant biomass of the undergrowth is important from a 
nutritional viewpoint for reclaimed post-mining site stand stabil¬ 
ity (Hiittl and Weber, 2001, Rnoche et al. 2002; Pietrzykowski 
and Socha 2011). 

Tree stand features, such as biomass and growth parameters 
(mean diameter and height), on the investigated sites may be 
more significantly altered by age and forest management opera¬ 
tions such as thinning directly affecting the density and competi¬ 
tiveness of trees. Thus, as mentioned in method description (see 
Section 2), the correlations with these parameters were not con¬ 
sidered. The values of MSQI correlated significantly with 
aboveground biomass of forest floor species. What follows is 
that MSQI well describes the growth conditions of vegetation 
that appears by way of natural succession under the canopy of 
the introduced tree stands. 

Proposal of final habitat classification and recommendation for 
species composition 

Based on the results of Tukey's HSD test, homogeneous subsets 
in terns of index values MSQI were distinguished (Table 6 and 
7). Mean MSQI values for three homogenous subsets (with low, 


medium, and high MSQI values) were used to determine the 
range of MSQI values for the grid with site classification (Table 
7 and Fig. 2). Based on this grid, the investigated soils develop¬ 
ing on different types of rock overburden may be grouped from 
the poorest to the most fertile soils. Thus, the poorest habitats 
were found on the Szczakowa sand pit on quaternary sands 
(SCZ-QS) (MSQI = 0.270), whereas the most fertile ones were 
found on Piaseczno spoil heap on quaternary sandy loams mixed 
with neogene clays (PIAS-QS+NC) (MSQI = 0.720) (Fig. 2). 

Eventually, the following predicted habitats were classified for 
the MSQI value ranges: below 0.403 - Coniferous Forests (CF); 
from 0.404 to 0.505 - Mixed Coniferous Forests (MCF); from 
0.506 to 0.602 - Mixed Deciduous Forest (MDF); above 0.602 - 
Deciduous Forests (DF). It must be noted however that the 
MSQI ranges for some soils overlapped within the groups. This 
was due to significant micro-habitat variability and mosaic-like 
character of sediments that form post-mining sites. In such cases, 
transition units (e.g., CF/MCF) may additionally be distin¬ 
guished. 

For individual groups of habitat classes distinguished accord¬ 
ing to similarity in substrates (parent rock), optimal species com¬ 
position for afforestation and conversion of pine monocultures 
would be recommended. 

For the poorest habitats of coniferous forest (CF) on quater¬ 
nary sands (e.g., Szczakowa sand pit SCZ-QS) and carboniferous 
deposits on unfertilized fragments of spoil heaps (Smolnica 
SMOL-NCF), the Scots pine should be still the main afforesta¬ 
tion species in the next generation of tree stands. For mixed co¬ 
niferous forest sites (MCF) on quaternary loamy sands (on Pi¬ 
aseczno spoil heap PIAS-QLS) and neogene sands (on Bel- 
chatow spoil heap BEL-NS) the dominant species should again 
be the Scots pine, but the addition deciduous trees (e.g., common 
birch and some sessile oak) is recommended. For transitional 
habitats, from mixed coniferous to mixed deciduous forest sites 
(MCF/MDF) classified on quaternary loamy sands (on Szcza¬ 
kowa sand quarry SZC-QLS), the dominant species in the next 
generation should be the Scots pine, but the introduction of more 
deciduous trees (e.g. sessile oak) is strongly recommended. 

Mixed deciduous forest sites classified on carboniferous de¬ 
posits of fertilized fragments of a spoil heap (Smolnica SMOL- 
CF), the conversion of pine monoculture is recommended. The 
main species on these sites should be deciduous trees (sessile oak 
with some other deciduous species, e.g., hornbeam, linden, 
common maple) and only a small percentage of the Scots pine. 
Transitional habitats, from mixed deciduous to deciduous forest 
sites (MDF/DF) were classified on a mix of quaternary sands 
with neogene clays (on Piaseczno sulphur mine spoil heap PIAS- 
QS+NC). The main species on these sites should be deciduous 
trees as well (e.g., sessile oak with admixture of hornbeam, lin¬ 
den, common maple) and a very small percentage of the Scots 
pine. 
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SCZ-QS SMOL-NCF PIAS-QLS BEL-NS SZC-QLS SMOL-CF PIAS-QS+NC BEL-QL 

Mine site and substrate variant 


Fig. 2. A grid of predicted groups of forest habitat units for the study sites and substrate (parent rock) variants, based on the range of MSQI: CF - poorest 
habitats of coniferous forests; MCF - mixed coniferous forest; MDF - mixed deciduous forest; DF - deciduous forest. 

Table 6. The result of Tukey’s HSD (honestly significant difference) test on the significance of mean MSQI (Mine Soil Quality Index) values for soil 
substrate variants 


Mine sites and substrate 

BEL-QL 

BEL-NS 

SMOL-CF 

SMOL-NCF 

SCZ-QLS 

SCZ-QS 

PIAS-QS+NC 

PIAS-QLS 

variant symbol 

M =0.625 

M =0.455 

M =0.573 

M =0.380 

M =0.486 

M =0.346 

M =0.608 

M =0.429 

BEL-QL 


0.003730 

0.861012 

0.000168 

0.025044 

0.000149 

0.999761 

0.000793 

BEL-NS 

0.003730 1 


0.083236 

0.526889 

0.990454 

0.131052 

0.010948 

0.996684 

SMOL-CF 

0.861012 

0.083236 


0.000974 

0.357413 

0.000229 

0.981651 

0.018603 

SMOL-NCF 

0.000168 

0.526889 

0.000974 


0.148320 

0.985062 

0.000222 

0.899769 

SCZ-QLS 

0.025044 

0.990454 

0.357413 

0.148320 


0.023257 

0.067971 

0.798833 

SCZ-QS 

0.000149 

0.131052 

0.000229 

0.985062 

0.023257 


0.000154 

0.410692 

PIAS-QS+NC 

0.999761 

0.010948 

0.981651 

0.000222 

0.067971 

0.000154 


0.002193 

PIAS-QLS 

0.000793 

0.996684 

0.018603 

0.899769 

0.798833 

0.410692 

0.002193 



1 0.007494 - marked differences by bold lettering are significant at p =0.05; M - mean of MSQI for mine soil and substrate variant. 


Table 7. Homogeneous subsets of soil-substrate variants in terms of MSQI values based on Tukey’s HSD test. 


Mine site and substrate vari- 

MSQI 



Homogeneous subsets 



ants 

(Mean for variant) 

1 

2 

3 

4 

5 

SCZ-OS 

0.346 






SMOL-NCF 

0.380 






PIAS-QLS 

0.429 






BEL-NS 

0.455 






SCZ-QLS 

0.486 






SMOL-CF 

0.573 






PIAS-QS+NC 

0.608 






BEL-QL 

0.625 






mean MSQI for homogenous subsets 

0.403 1 

0.438 

0.505 

0.555 

0.602 


1 0.403- marked mean values by bold lettering for homogeneous were subsets used as the ranges for predicted units in forest habitat classification 


^ Springer 





















Journal of Forestry Research (2014) 25(1): 63-74 


73 


The most fertile habitats of deciduous forest sites (DF) were 
classified mostly on Quaternary loam (on ‘Belchatow’ spoil heap 
BEL-QL) (Table 6). On these deposits the conversion of pine 
monoculture and afforestation with hard wood species is strongly 
recommended. 

Conclusions 

Reclaimed mining soils developing on post-mining sites are 
made up of very diverse lithological deposits that exhibit consid¬ 
erable variability. The main factor that affected the variability of 
MSQI for quaternary and neogene soils (BEL-QL, BEL-NS, 
SCZ-QLS and SCZ-QS PIAS-QS+NC and PIAS-QLS) was soil 
texture, which resulted in the value of soil texture index (1st). 
Moreover, in the case of soils developed on shales and Carbonif¬ 
erous deposits of Smolnica spoil heap (SMOL-CF and SMOL- 
NCF), the features that also had an impact were fertility and the 
resulting sub-indices, including nutrient availability index, Ina, 
and phosphorus availability index IP. 

On this basis, it may be concluded that, with little lithological 
diversification of sediments on the Smolnica spoil heap, the im¬ 
pact of mineral fertilization was highlighted. The analysis and 
habitat classification results lead to the conclusion that prediction 
of plant habitat based solely on lithology and genesis of sedi¬ 
ments made of parent rock, in the case of developing mine soils, 
provides results that are too general. Often the predicted potential 
fertility of weathering rock overburden based on an analogy with 
the parent rocks of natural soils, (e.g., shale, and sandstone) may 
be erroneously interpreted and ultimately give overestimated 
final results of habitat classification. The conducted habitat clas¬ 
sification, especially the examples provided, show that the de¬ 
veloped MSQI index is universal and can objectively assess the 
quality of mine soils in contrast to assessments based solely on 
soil properties that are typically highly variable. 

Moreover, the reliability of the developed MSQI index is con¬ 
firmed by statistically significant correlations with the tested 
features of plant communities from succession (ecological indi¬ 
cator values and plant biomass). This indicates the accuracy of 
the selected components of soil quality assessment included in 
MSQI and determines fertility of the soils and ability to meet 
nutritional requirements of plants. As a result of the habitat clas¬ 
sification and an indication of the trend that habitats will follow 
in the process of becoming similar to natural forest habitats, the 
MSQI index can be useful in designing species composition and 
pine monoculture transfonnation in the next generation of “new 
forests” on mining sites. The MSQI can by applied to describe 
the variability of mine soils and to classify habitats developed on 
mine overburden deposits, similar to the investigated substrates, 
which are dominant in this part of Europe. 
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